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Abstract 

We study the transition of a granular packing from liquid to solid 
bonding in the course of drying. The particles are initially wetted by 
a liquid brine and the cohesion of the packing is ensured by capillary 
forces, but the crystallization of the solute transforms the liquid bonds into 
partially cemented bonds. This transition is evidenced experimentally by 
measuring the compressive strength of the samples at regular intervals of 
times. Our experimental data reveal three regimes: 1) Up to a critical 
degree of saturation, no solid bonds are formed and the cohesion remains 
practically constant; 2) The onset of cementation occurs at the surface and 
a front spreads towards the center of the sample with a nonlinear increase 
of the cohesion; 3) All bonds are partially cemented when the cementation 
front reaches the center of the sample, but the cohesion increases rapidly 
due to the consolidation of cemented bonds. We introduce a model based 
on a parametric cohesion law at the bonds and a bond crystallization 
parameter. This model predicts correctly the phase transition and the 
relation between microscopic and macroscopic cohesion, 
keywords: Unsaturated granular media, capillarity, cementation, DEM, 
compressive strength 
PACS: 45.40.-f; 45.70.Mg; 81.05. Rm 

1 Introduction 

Granular materials occur very often in nature and industrial applications with 
liquid or solid bonds between the particles. The effect of capillary bonding on the 
mechanical behavior of granular materials is of primary importance in powder 
technology [1] [2j [3] and transformations of geomaterials [U |6]. Capillary 
bonding has been extensively studied in the past, and several models of capillary 
cohesion have been proposed [7J [HJ [S] . The scale- up of capillary forces to the 
macroscopic cohesion of granular media has been investigated by means of both 
experiments and discrete element numerical simulations [131 1141 115] . 
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In the same way, solid bonding is found in various materials such as sedi- 
mentary rocks (sandstones, conglomerates and breccia) [24] , biomaterials such as 
wheat endosperm (starch granules forming a compact structure bound together 
by a protein matrix) [551 HH], and geomaterials like mortars, concrete and as- 
phalt (aggregates of various sizes glued to each other by a cement paste) [2"T] . 
The cohesive behavior due to solid bonds between the particles and the rup- 
ture of cemented granular materials have been recently investigated by several 
authors [HUTU [25]. 

There is a broad class of processes where the nature of cohesive bonds evolves 
from solid to liquid and vice versa. The liquid is generally water with various 
types of impurities and / or dissolved minerals. As a result of reactive transfers or 
phase transitions, the liquid may influence the mechanical properties by dissolu- 
tion of cemented joints or crystallization of solutes. The mechanical properties 
can thus evolve in response to external hygrothermic conditions [8]. In phar- 
maceutical industry, for example, the process of wet granulation is widely used 
to manufacture tablets. A liquid composed of solvent(s) and solute(s) is mixed 
with dry powder in order to obtain a homogeneous mixture that would segre- 
gate in the absence of the liquid. The formation of aggregates in the process of 
mixing has been a subject of extensive investigation Q] [H [3] . The solid tablets 
are obtained by drying: the solvent evaporates and solid bridges are formed 
between grains by crystallization, thus ensuring the macroscopic cohesion of the 
tablets. 

In this paper, we are interested in the evolution of the macroscopic cohe- 
sion of a wet granular material during the evaporation of a solvent (water) and 
thus crystallization of the solute (sodium chloride). At low liquid contents, the 
water is in the form of capillary bridges between particles. The hydrous field 
corresponding to this range of liquid contents is often described as "pendular" . 
The attractive capillary forces in the presence of liquid bridges between parti- 
cles endows the material with macroscopic cohesion [8[ [29] . Under the action of 
the surrounding hygrothermic conditions, the water evaporates and the solutes 
crystallize to form solid bonds between grains. The cohesion is initially of cap- 
illary origin but begins to increase as soon as the first solid bonds are formed at 
the free boundary of the sample and a solidification front propagates from the 
boundary to the center of the sample. 

We measure the compressive strength of several samples of glass beads and 
sand, prepared in the pendular state by mixing with a saturated brine, in the 
course of drying. Technically, it is not trivial to follow simultaneously the evolu- 
tion of the internal state due to the evaporation of the brine and the evolution of 
its mechanical behavior. Moreover, the two limits of purely liquid bonding and 
purely solid bonding have different behaviors and require specific experimental 
tools and precautions. For these reasons, the experimental campaign involves 
the preparation of many reproducible samples subjected to identical drying con- 
ditions, each sample being destroyed when subjected to simple compression test. 

We show that the evolution of the cohesion is governed by three time scales 
corresponding to the evaporation of a single capillary bridge, the propagation of 
the cementation front from the surface to the center and the aging of cemented 
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bonds. These time scales are identified from the experimental campaign together 
with a theoretical model based on a a simple representation of the evolution of 
the debonding force between the particles. This model fits well the compressive 
strength as a function of an order parameter representing the state of crystal- 
lization. 

In the following, we first describe in section [5] the experimental setup, the 
procedures of sample preparation and the main results regarding the evolution 
of the samples and mechanical strength in the course of drying. In section [31 
we present our model and compare its predictions with the experimental data. 
We conclude with a summary and scopes of this work. 

2 Experiments 

2.1 Experimental setup 

The experimental tests were designed to measure the compressive strength of 
granular samples during the process of the drying and deposition of a solute 
at the liquid bridges between the grains. Two types of materials were used in 
these tests: (1) a model material composed of glass beads and (2) a natural 
material: Ventoux sand. In both cases, the grains were first washed and dried. 
Then, they were sieved to keep only the grains with diameters in the range 
from 0.4 mm to 0.8 mm. The dry grains were homogeneously mixed with a 
small amount of saturated brine containing 35.6 g of NaCl for 100 g of water in 
controlled environment (temperature, humidity). This maximum concentration 
varies with temperature jTTJ [TH] . 

The samples (diameter 25 mm and height 17 mm) were prepared inside a 
cylindrical mould. At this stage, the samples are particularly delicate to han- 
dle since their mechanical strength is solely ensured by weak capillary forces 
acting between the grains. Because of this brittleness, particular experimen- 
tal precautions were necessary in order to avoid damage to the samples upon 
unmoulding. 

Each sample is initially a three-phase material including 

• a solid phase (denoted by S) composed of grains (glass beads or sand) 
with a total mass rrig rains , 

• a liquid phase (denoted by L) composed of water of mass mfj 2 o an d 
dissolved NaCl of mass rn^ aCl . 

• a gaseous phase composed of air and water vapor with negligible mass 
compared to the two other phases but which carries the kinetics of drying. 

Thus, the total mass of the sample at this initial state is to^ = ™^ 2 o ~^ m NaCi + 

777 & 

grains ' 

The water evaporates at the surface and causes thermodynamic disequilib- 
rium of the solution. The equilibrium is recovered through the crystallization of 
NaCl and thus the formation of solid NaCl deposits of a total mass rnff aCl . A 
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Figure 1: (a) Photograph of the sample in the process of drying; (b) The press 
used for the unconfined compression tests. 



crystallization front propagates thus from the boundary towards the center of 
the sample. The evolution of a sample is expected to be governed by the ratio 



Ic, = 



m NaCl 



m NaCl 



m NaCl 



(1) 



to which we refer as the global crystallization index of the sample. 

Three series of samples were prepared with three different values of the initial 



liquid content Wn= {mj 



,)/m rams : 3 %, 5 % and 7%. These values 



l H20 + m NaCl ) I m gram > ' 

were selected in the hydrous field corresponding to the "pendular state" where 
liquid bridges between grains ensure the mechanical integrity of the sample. 
The samples were weighed just after preparation, and left to dry at constant 
temperature T = 20° C and relative humidity RH = 43 % for drying times 
varying from 15 minutes to 18 hours; see Fig. QJa). 

At regular intervals of time, the samples were weighed (mass to/), and then 
subjected to unconfined vertical compression up to failure. The tests were car- 
ried out by means of a "low capacity" press, which allows us to test samples 
with low dimensions and forces up to 50 N with an accuracy of 0.01 N; see Fig. 
|ljb) . The difference mj-m t yields the mass of the evaporated water and thus 
that of the crystallized mass fraction of NaCl. Since all samples are prepared 
under similar conditions, they are assumed to follow the same evolution in time. 
This allows us to reconstruct the kinetics of crystallization. 



2.2 Experimental results 

The evolution of the compressive stress <j y at yield is shown in Fig. [2] as a 
function of the global crystallization index I c for glass beads and sand, as well 
as from a theoretical model described in section [3] All the experimental data 
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experiments: 




Figure 2: Compressive yield stress a Y , normalized by the yield stress a\ due 
only to capillary forces (yield stress in the first regime) as a function of the 
global crystallization index I c in experiments. The solid line corresponding to 
theoretical model (see section 3) is plotted as a guide to eyes. 



collapse nearly on the same curve for glass beads and sand. We distinguish three 
regimes in the evolution of the cohesion. As long as I c is below a threshold 
It, the yield stress is independent of I c . In this regime, the amount of the 
crystallized salt at the boundary layer of the sample, where the whole process 
is initiated due to exchange with the atmosphere, is insufficient to allow for the 
formation of solid bonds between the particles. Therefore, the yield stress a Y 
in this regime is only due to the capillary bonds between the particles. In Fig. 
[H we see that a Y is independent both of the initial liquid content wn and of 
the type of grains. In fact, both experiments and theoretical analysis suggest 
that in a homogeneously wetted granular material in the pendular state, the 
capillary cohesion is practically independent of the liquid content but crucially 
depends on the number density of the capillary bonds [29] . 

The yield stress increases from I c ~ 0.9 due to the partial cementation of the 
first bonds at the boundary of the sample and the propagation of a solidification 
front towards the center of the sample. The yield stress increases rapidly in this 
regime. A photograph of the cementation front in a partially dried sample is 
displayed in Fig. [3] The bonds in the central part of the sample were of capillary 
type and thus the grains were easily removed. The remaining grains form a solid 
crown due to cemented bridges between the particles. Hence, the thickness of 
this crown represents the distance run by the solidification front. It should be 
noted that the solid bonds behind this front are only partially cemented. At 
this stage, each solid bond is composed of a crystallized superficial layer with 
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Figure 3: Photograph of a partially dried sample in the second regime. 




Figure 4: Schematic representation of the evolution of bridge: from a liquid to 
a solid bridge. 

liquid brine inside. As drying continues, the cemented bond gets thicker and 
stronger until the whole bond solidifies. This process is illustrated in Fig. |U 

All bonds are cemented when the center of the sample is reached by the 
cementation front. In this third regime, a Y rises linearly to a level as high as 35 
times the capillary yield stress due to the consolidation of the cemented bonds 
throughout the system as a result of continued crystallization of the brine at each 
bond. These high values of compressive strength as I c — > 1 has two different 
origins. On one hand, the debonding force of a solid NaCl bond is nearly 7 
times that of a capillary bond. On the other hand, unlike a liquid bond which 
involves only a tensile strength in addition to sliding friction, a solid bond has a 
rolling (or bending) strength as well as a tensile strength. This rolling strength 
is crucial for the compressive strength of cohesive granular materials [3TJ [23] . 

Figure [5] shows the evolution of I c with time t. This evolution reflects the 
kinetics of drying with the transport of water in the gas phase. I c begins to 
increase as soon as the sample is exposed, but the first partially cemented bonds 
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Figure 5: Evolution of the global crystallization index with time; solid line: 
experimental results, dotted line: model equation ([2]). The error bars represent 
the standard deviation for five independent experiments. 

are completely formed only at I c = It — 0.9. Then, I c increases mainly due to 
the formation of the first (partially) cemented bonds followed by the propagation 
of the cementation front into the sample. When this front reaches the center of 
the sample, I c increases only due to the aging and consolidation of the cemented 
bonds. Although the kinetics of drying is not exactly the same in each regime, 
the time evolution of I c is quite well fit by an exponentially increasing function: 



where td is a characteristic drying time. This form suggests that the drying rate 
is proportional to the residual liquid. 



The evolution of the compressive strength may be understood on quantitative 
grounds by considering 1) a local cohesion law involving the phase change from 
liquid to solid at the level of a single bond, 2) the relation between the bond 
strength and the macroscopic strength, and 3) the kinetics of drying both at 
the surface and in the bulk. In the following, we present a model that is based 
on these ingredients and correctly captures the macroscopic evolution as shown 
by the fitting form in Fig. [2] 



h{t) = 1 — e 



(2) 



3 Model 
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Figure 6: The parametric adhesion law. 



3.1 Single-bond behavior 

The phase change of a bond from liquid to solid is controlled by the amount of 
the deposited solute. Hence, the natural control parameter for phase transition 
at this scale is the bond crystallization index i c defined by 

• _ m NaCi( hcmd ) 

m NaCl ( b0nd ) + m NaCl (bond) 

This index depends on the history of the bond and is thus a function i c (r, t) of 
the radial position r of the bond and the time t elapsed since the beginning of 
drying. We assume translational invariance along the axis of the sample. The 
index i c varies between for a liquid bond and 1 for a fully cemented bond. The 
transition from liquid to solid takes place at a particular value i c = it which 
corresponds to the deposition of the first percolating shell of crystallized solute 
across the gap or at contact between two particles. In this transition, the tensile 
force threshold f Y (the debonding force) increases from that of a liquid bond 
fY to that of a fully cemented bond fg . We have f Y = fY for < i c < it and 
/ = fg for i c = 1. Assuming a linear evolution between these two limits, we 
get the following parametric adhesion law: 

f = fl + H(i c - i T ) ^(/J - fl) (4) 
i — it 

where H is the Heaviside step function. A plot of this relation is shown in Fig. 

m 

Remark that f Y in equation 2] is the force threshold and the actual value of 
the normal force /„ depends on the force law expressing a force-displacement 
relation. This relation is different for a liquid bond and a partially or fully 
cemented bond. Several expressions have been proposed for the capillary force 
fl as a function of the gap S, bond liquid volume Vb and surface tension 7 
[SI HH] ■ The capillary force has its largest value at the contact between two 
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particles and declines exponentially with the gap: 

f L = 2TT 1C os9 P e- 5 /\ (5) 

where A is a characteristic length depending on V& and the particle radii, p 
is a reduced length depending on the particle radii and 9 is the contact angle 
[HJ [32 ■ A capillary bridge is stable as long as S < S max given by [33] 

5 max = (i + ^K 73 - (6) 

This distance is small (compared to A) in the pendular state and in practice the 
debonding force can be approximated by 

/r = 27r 7 cos0p (7) 

This force is independent of the bond liquid volume VJ, . 

The fully cemented bonds have a brittle behavior with a debonding threshold 
fg ~ 7fY for the three values of wn used in our experiments. The behavior 
is also brittle for the partially cemented bonds with a threshold which is as- 
sumed here to vary linearly with the bond crystallization index according to 
equation [4] Obviously, the cemented bonds exhibit also a shear strength and a 
rolling strength in addition to the tensile strength considered above. This fea- 
ture, which differentiates cemented bonds from liquid bonds, is not considered 
explicitly in the single-bond behavior but will be taken into account through a 
Coulomb cohesion parameter c at the macroscopic scale. 



3.2 From single bond to RVE 

A schematic representation of our model is shown in Fig. [7j It consists of two 
scales: 1) The microscopic or local scale is defined by the parametric adhesion 
law f Y {i c (r, t)}; 2) The macroscopic or global scale is characterized by the 
relation <r Y {I c (t)}. Given the local adhesion law, the macroscopic behavior 
can be determined from the relationships between i c and I Cl on one hand, and 
between f Y and er y , on the other hand. We first consider a Representative 
Volume Element (RVE). Then, solve the problem in the presence of hydric 
gradients in the drying process by introducing a simple propagation equation 
for the local crystallization index. 

We assume a homogeneous distribution of the liquid bonds in the initial 
preparation of the sample. We also assume that the total mass of NaCl is the 
same for all bonds. Then, the macroscopic crystallization index is simply the 
volume average of the bond crystallization index: 

R 

hit) = (*c(r,t))v = y / 2nri c {r,t) dr (8) 
o 

where h and R are the height and radius of the sample, respectively, and V — 
irR 2 h is its volume. According to equation for a RVE and in the absence 
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IS) ° Y (t) 
i c (r,t) * > f Y (r,t) 

Figure 7: Diagram of microscopic model and upscaling method. 



of hydric gradients, we have It — ir since all bonds reach simultaneously this 
point. During drying, the cementation is initiated at the surface of the sample 
for i c = %t and we have i c < it everywhere in the bulk. 

In order to propose a relation between f Y and a Y , we need to postulate the 
deformation mechanism, the macroscopic failure mode and the general expres- 
sion of the stress tensor in a granular material as a function of the bond forces. 
The stress tensor cry is the first moment of the bond forces fi averaged over the 
control volume V [23l [30] : 

c»j = n b (fdj)v (9) 

where n b is the number density of the bonds and lj is the j-component of the 
branch vector joining the centers of the two particles at contact. The averaging 
runs over all contacts belonging to the control volume. Assuming that the 
branch vector lengths and forces are not correlated and for a medium with 
uniformly distributed branch vectors in all directions, the normal stress a n is 
given by 

°n = \n b (l)(f n ) (10) 

where f n is the normal bond force. We have (ft) ~ as a result of the balance 
of force moments on each particle. 

According to equation [TU1 the maximum tensile force a™ ax in a RVE is 
obtained by replacing /„ by f Y : 

< MX = \n b (t)(f Y ) (11) 

Assuming that the Mohr- Coulomb failure criterion holds in the tensile regime, 
the Coulomb cohesion c is thus simply given by 

c = K WI -^ b tan^)(/ y ) (12) 

where fi is the friction coefficient and tp is the friction angle. 

Finally, we need the relation between the axial stress a\ and the cohesion 
c in order to evaluate the compressive strength. Under the conditions of axial 
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symmetry, we have a% = 03. Let p = (<7i + 203) /3 and g = (<7i — 0-3 )/3 be the 
mean stress and stress deviator, respectively. At failure, we have 



- = sin p H — cos (13) 

p 3 — sin tp \ p J 

In our experiments the confining stress is zero: a% = (73 = 0. Hence, from 
equation (fl"3"| . at failure we have 

CTi = a 5 = 2c 1 — (14) 

1 — sin 95 

Finally, the substitution of c from equation 1121 yields 

3 1 — sin <p 

Equation (1151) shows that the compressive strength is proportional to the 
microscopic force threshold and to the number density of the bonds. Notice 
that rib is simply half the average number of bonds per particle divided by the 
free volume, i.e. the mean volume V p of a Voronoi cell surrounding the particle: 

* - k (16 > 

By definition, given the solid fraction ip of the sample, V p (f> is the mean volume 
of a particle. Thus, 

where d is the particle diameter. For nearly spherical particles, we may set 
(£) ~ (d). From equations (fT5|) . (fTl?)) and (IT71) . we get the following equation in 
which the prefactor to f Y is expressed in terms of the familiar parameters: 



Y 



a = —z 



2 sin ip (d) Y 



7r 1 — sin (f (d 3 ) 



(f ) (18) 



For the typical values z — 6, ip — ir/6 and 4> = 0.6 and assuming monodisperse 
particles, we get a Y ~ 2f Y /d 2 , which is a simple relation. 

3.3 Compressive strength 

From the relations jSJ) and ([TBI) together with the parametric adhesion law (|4|, 
it is straightforward to predict the macroscopic behavior. We first take the 
average of both sides of equation : 

(f Y ) = fl + 4 ~ Sl - fr) g (*° - *r)) (19) 
1 — 
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where the averaging runs over the whole sample. Inserting this expression of 
the mean force threshold in equation (fT8|) , we get 



a = —z 



{ft + f ~T^ {{lc ~ 1t)h (lc - iT)) } (20) 



7r 1 — sin y> {d 3 } 

The average ((i c — ir)H (i c — ir))v in the right hand side of equation (|20p 
depends on the function i c (r, i) and will be evaluated below in the presence of a 
hydric gradient. Here, we consider the first and third regimes identified in the 
experiments (see section |2~2|) . The first regime corresponds to the initial stage 
of drying where i c (r,t) < it for all bonds and no cemented bond is formed. 
Hence, according to (|20|) . the compressive strength in this regime is given by 

- y =*r=-*V^Jr/z (21) 

7r 1 — sin ip \d 6 ) 

We see that o Y in this regime is constant and independent of I c . The purely 
liquid-bond strength a Y can be used to scale a Y in equation (|20p: 

° =1+ JsHl_ {{1c _ 1t)h{1c _ 1t)) (22 ) 

In the third regime, the cementation front has reached the center of the 
sample and we have it < ic(f,t), and thus H(i c — it) — 1 everywhere. As a 
result, from equations © and (|2"2"j) we get 



Y fY I fY i 

^^1+ /s / /l - (Jc-vr) (23) 



This equation shows that in the last stage of drying, corresponding to the con- 
solidation of all bonds, the compressive strength increases linearly with I c . This 

linear behavior and the predicted value of the coefficient ? {_f T — — 30 for 
fs/f Y — 7 and %t — 0.9 are consistent with the data shown in Fig. [2] within 
our experimental precision. 



3.4 Cementation front 

The second regime is a transient beginning with the cementation of the bonds 
located at the surface and ending when all bonds are partially cemented. In this 
regime, the bonds at a radial position r < rx(t) are liquid whereas the bonds 
with rx(t) < r are partially cemented. The cementation front rx(t) decreases 
from R at the end of the first regime to at the beginning of the third regime. 
By definition, we have i c (r = Vt, t = iy) = it- 

For the evaluation of the compressive strength in the second regime, we need 
the function i c (r, t). By definition, 1— i c (r, t) is the proportion of the liquid brine 
in the bond, which is conserved in the bulk: The time variation dt{l — i c (r, t)} 
of the liquid in the bonds located at r is compensated by a flux J in the vapor 
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phase. Assuming that J is proportional to the gradient d r {l — i c (r, t)}, we arrive 
at a diffusion equation for i c : 



dt \ dr 2 r dr _ 

where invariance by rotation around the axis and by translation along the axis 
are assumed, and D is a diffusion coefficient. 

Evaporation takes place only radially and at the free surface r = R. The 
exponential increase of I c observed in the experiments (see Fig. [5]) suggests a 
simple exponential time evolution of i c {r — R,t): 

i c (R,t) = 1 - e-W^ (25) 

The characteristic time represents the time necessary for the onset of cemen- 
tation in a bond when i c = vr. We have ix = 1 — . Hence, equation (|25p 
can be rewritten in the following form 

i c (R,t) =1-(1- i T ) t/tT (26) 



Equation (f24[) can be solved with (|26j) as its boundary value for an initially 
homogeneous system in which i c (r, t = 0) = for all r ^ R. The solution is [32] 

^ *) = {i - a - ^ > { i - 1 E e " p t^y } (27) 

( a n Ji(a n R) J 

where Jo and J\ are Bessel functions of the first kind of order zero and one, 
respectively, and a n is the rfi 1 root of the equation aJ\(aR) — Jo(aR). The 
solution (eq. [57|) is displayed on Fig. [5] 

From equation ([27]) we can calculate the global crystallization index defined 
by©: 

9 7 { 4 °° p - Dq 4* 1 

IS) = - 2 J ri c (r,t) dr = {1 - (1 - fa.)*/**} 1 - ^ £ — — (28) 



o 



n=l 



The macroscopic transition index It is the value of I c at t = tr- Hence 

4 

R 2 



n=l n ) 



This relation implies that It < *t- Given the high value of ty, in practice we 
have It — It- 

Equation (1271) is an implicit equation for the cementation front rxif) defined 
by the condition i c (rr,i) = it- The knowledge of rx{t) allows us to determine 
the total compressive strength a Y as the weighted mean of u\ in the region 
r < tt and erg in the region tt < r: 
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Figure 8: Evolution of the local crystallisation index i c in the bulk with time. 
Plugging the expressions (|2"TT) and (|2"2"|) into equation ([50)1 yields 

a? 1 + i 1 - "R^J l-vr {h ~ lT) (31) 

where it may be replaced by It from equation 1291 

Equation (J2J) holds only in the second regime, but we retrieve the expression 
of Og in equation (j2"2")l for the third regime by setting r-r = as well as the first 
regime by setting tt — R (equation (|2"2"10 . The nonlinear increase of a Y in the 
second regime with I c is given by the factor 1 — r T /R 2 . This corresponds to a 
first-order phase transition at the macroscopic scale for I c as control parameter 
and <j y as state parameter. 

The model described in this section is able to predict the global crystal- 
lization index I c {t) as a function of time (equation (|T)>) and the compressive 
strength cr Y (J c ) as a function of I c (equations (|5T|) . (j2"2"j) and (|2"2"1) '). The plots of 
these functions for the experimental parameters are displayed in figures [5] and 
[21 The predicted behavior is thus fairly good within the available experimental 
precision. Apart from the radius R of the sample, for these fits three parameters 
were used from experimental measurements, namely £x> D and iy. The good 
estimation of a\ provides a strong support for our model. 

4 Conclusion 

In this paper, we introduced a simple experimental setup allowing us to ana- 
lyze the transition from capillary cohesion to solid cohesion induced by surface 
drying and the progressive crystallization of a liquid brine wetting initially the 
particles. This transition was evidenced experimentally by measuring the com- 
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Figure 9: Influence of the local transition index on the compression strength as 
a function of the macroscopic crystallisation index I c . 

pressive strength of the samples in the course of time, and modeled by intro- 
ducing the relevant variables both at the scale of a bond and at the macroscopic 
scale. Three regimes were evidenced: 1) the capillary regime as long as the 
degree of saturation remains below a critical value; 2) the transition regime, 
corresponding to the initiation and propagation of a cementation front spread- 
ing towards the center of the sample, and 3) the cemented regime, where all 
bonds are partially solidified but continue to strengthen due to further crystal- 
lization. The predicted behavior by our theoretical model is in good agreement 
with the experimental data. This model, based on an evolutive cohesion law 
at the bonds, takes into account the kinetics of crystallization and the granu- 
lar texture. It correctly predicts the compressive strength in the capillary and 
cemented regimes. 

The framework of the model presented in this work is quite general. It 
can be applied to granular materials subjected to evolutive force laws or phase 
compositions. For example, the degradation of a geomaterial under the action 
of progressive dissolution of solid bonds in a humid environment bears strong 
similarity with the problem treated in this paper. 
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